Spin transport properties of a quantum dot coupled to ferromagnetic leads with 

noncollinear magnetizations 



o 
o 

(N 



oo 



Hao Zhang^, Guang-Ming Zhang^, and Lu Yu'^'^ 
^Institute of Theoretical Physics, Chinese Academy of Sciences, 100190 Beijing, China; 
^Department of Physics, Tsinghua University, 100084 Beijing, China; 
''institute of Physics, Chinese Academy of Sciences, 100190 Beijing, China. 
(Dated: March 18, 2009) 

A correct general formula for the spin current through an interacting quantum dot coupled to 
ferromagnetic leads with magnetization at an arbitrary angle 6 is derived within the framework of 
the Keldysh formalism. Under asymmetric conditions, the spin current component Jz may change 
sign for < S < vr. It is shown that the spin current and spin tunneling magnetoresistance 
exhibit different angle dependence in the free and Coulomb blockade regimes. In the latter case, the 
competition of the spin precession and the spin-valve effect could lead to an anomaly in the angle 
dependence of the spin current. 
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I. INTRODUCTION 

In the new field of spintronica^, the magnetic proper- 
ties of quantum devices control the transport properties 
via the electron spin, for example, the tunnel magnetore- 
sistance (TMR) in ferromagnetic tunnel junctions. The 
high magnetoresistance in a TMR device is due to the 
spin- valve effect, namely, the resistance strongly depends 
on whether the magnetization of the two ferromagnetic 
electrodes are parallel or antiparallel. By switching the 
magnetization of one electrode with respect to the other, 
the charge current is modulated by the relative angle of 
the two magnetic moments. With the magnetic tunnel- 
ing injection technique, a pure spin current can be gen- 
erated and detected experimentallyi^i^ This substantial 
progress in experiment makes it feasible to investigate 
the spin transport properties in mesoscopic systems. 

To study the spin-dependent transport properties, a 
device setup of a quantum dot (QD) coupled to ferro- 
magnetic leads has been proposed."* In such a geome- 
try, the charge current can be spin polarized and can 
induce a pure spin current. However, up to now, most 
of the previous works were devoted to the charge trans- 
port properties, not to the study of the spin current 
itself. Moreover, the main focus was on the charge 
transport on a QD coupled to two ferromagnetic leads 
with coUinear magnetizations^ii^i^iiSiiiii^ii^iiiii^ii^ 
while less attention was given to the noncollinear 
alignment^OSis^^^^^^^^Braun et al^ gave 
an expression for the spin current through the left tunnel 
barrier. However, they did not derive an adequate uni- 
fied formula for the spin current through the two tunnel 
barriers and did not actually consider the spin current in 
a steady state. The ^-component of the spin current de- 
fined as a difference between the spin-up and spin-down 
contributions to the charge current was considered by Mu 
et alM for the noncollinear case. Unfortunately, these 
authors did not properly take into account the difference 
of the spin quantization axis for the two leads, so their 
result is correct only for the parallel case. 



Recently, Rudzihski et al^ studied the charge cur- 
rent through a quantum dot coupled to noncoUinearly 
polarized ferromagnetic leads. They found that the 
current- voltage curve reveals typical step- like characteris- 
tics. They also found that the spin precession is enhanced 
by the Coulomb correlations and strong spin polariza- 
tion of the leads. Moreover, the relationship between the 
charge current and the angle of the magnetization con- 
figurations of the electrodes has been studied by Zhou et 
al?^ These authors concluded that the angle dependence 
of the electric current in the free regime varies monotoni- 
cally from the parallel to antiparallel alignment, while in 
the Coulomb blockade regime it varies nonmonotonically. 
However, authors of both references did not consider the 
spin current in this general configuration. 




FIG. 1: Sketch of the system configuration. QD is coupled to 
two ferromagnetic leads with magnetizations Ml and Mu at 
an angle 9. 

In this paper, we first derive an exact general formula 
for the spin current through a QD coupled to noncollinear 
ferromagnetic leads, starting from the Heisenberg equa- 
tion for the spin operator in terms of the Keldysh Green's 
functions^^ (Sec. II). To the best of our knowledge, this 
general formula of the spin current for this class of devices 
is derived for the first time. It should play a similar role 
as its charge counterpart derived earlier^i^i^ Then, by 
using the equation-of-motion technique with the Hartree- 
Fock decoupling scheme, the spin current is obtained as 
a function of the bias voltage and the angle 6 of the mag- 
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netization configurations of the leads (Sec. III). Further- 
more, the spin current and the spin tunneUng magne- 
toresistance (STMR) are calculated numerically in both 
free and Coulomb blockade regimes (Sec. IV). The in- 
terplay of the spin precession enhanced by the Coulomb 
repulsion and the spin valve effect gives rise to anomalous 
behavior of the angular dependence of the spin current 
anticipated in the Coulomb blockade regime. 

II. GENERAL EXPRESSION FOR THE SPIN 
CURRENT 

The system considered in this paper is schematically 
shown in Fig.l, and it consists of a single-level quan- 



tum dot coupled to two ferromagnetic metallic leads by 
tunneling barriers. The magnetic moment M of the left 
electrode is pointing to the z-direction, while the moment 
of the right electrode is at an angle 6 to the z axis in the 
x-z plane. We will use the local and global quantization 
axes to describe the electron spin. The local quantiza- 
tion axes are determined by the local spin polarization in 
the leads, while the global axes are the local basis in the 
left electrode. The corresponding model Hamiltonian is 
given by^ 



h.c 



W,a;a=L,R 



E 



^k.if 4,+,fl cos - - cl _ jj sin - Ut + TkM 4 cos - + cj^ sin - Ui + h.c 



(1) 



where the spin projection on the local axes is denoted as a = ±, ek,CT,a = ek,a + uMa is the single-electron energy in 
the a-th electrode, and cj^ ^ ^ and Ck,a,a correspond to the creation and annihilation operators, respectively. Similarly, 
the spin projection on the global axes is denoted as 7 —H, dly and dj are the creation and annihilation operators of 
the electron on the quantum dot with energy Cd- 

For simplicity, we can rewrite the model Hamiltonian into a compact matrix form 



\i,a=L,R 



U 



Dtf))'- (Dta.f))' + (c^Tk,aRlD+/i.c.) , (2) 
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where we have introduced the Nambu spinors and two useful matrices 



C'k, + ,a 
C'k, — ,Q 



dl 



ek, + ,a 

ek,-,a 



(3) 



with Ol = for the left lead and Br ^ 9 for the right lead. When the spin operators of the two leads are considered 
^k 



Sa = {h/2) Q,o'aCk,Q, the spin matrices are given by 



1 
= i 1 



-i 

1 



sm f cos f 
cos — sin t 



V V z 

'<^R = <^L^^R 



1 

-1 y 

cos 9 — sm t 
— sin 9 — cos t 



(4) 



From the Heisenberg equation, we can calculate the spin current J a = (Jq) from the lead to the dot 
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h = J[S„,H] = ^^Tr(ct^^a„tk.aRiD-DtR„t*^^a„Ck,a 



(5) 



Moreover, by introducing the Keldysh Green's function matrices 



(<+,a(0)rfiW) {CL.jO)di{t)) 



,G<(t) = z 



(dj(o)dT(t)) (4(o)dTW) 

{d\iO)d^{t)} {d\iO)d^it)) 



(6) 
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we can further rewrite the expectation value of the spin current as 



du ^ 
2^ 



Re 



Tr(G< (k,c^)a„tk,oRi 



(7) 



where (k, w) is the Fourier transform of (k, <). Since the ferromagnetic leads are noninteracting, we obtain 
the Dyson equation for G^ (k, uj) in terms of the Green's function matrices for the local dot electrons, 



G< (k, u) = G^ (lo) R„n,«g< (Kto) + G< (lo) RaTl,aSl (k,c.) , 



(8) 



where 



ia (k,w) = 27ri/a (w) 
/ 1 



(5 (cj - ek, + ,a) 

^ 

1 



with fai(^) = [1 + exp((jj — ^a)/ (kBT)] ^, /iL = —eV/2 and /i^j = ey/2. Inserting these expressions into 
current formula, we obtain the spin current as follows: 



-Re Tr zr„ H 



2/„ (w) G^ {u) + G< (c^) - iP 



dE G< (ii;) 

TT E ^ LO 



the spin 



(9) 



where the integral is taken as the principle value and 







ek,iT,( 



(10) 



Since this system is quasi one-dimensional, different from the spin Hall systems^i^i^ in which the spin-orbit coupling 
is essential, we do not take into account those spin flip processes due to the spin-orbit coupling. So we do consider 
the spin current through QD as a continuous and conserved quantity. The steady state is realized in the system 
through the scattering process which is similar to the charge transport. As far as we understand, no one has studied 
the detailed relaxation process within the QD. In a steady state, the spin current is uniform, so J l = — J_r. Thus, we 
can symmetrize the spin current as J = (Jl — Jfl)/2 which is similar to the operation performed on the expression 
for the charge current^i^i^. The general expression for the spin current is then given by 



i/^^Re| 
2 7 27r 1 



Tr 



ilil 



Braun et al^ gave an expression for the spin current 
through the left tunnel barrier, but they did not derive an 
unified formula for the spin current through the left and 
right tunnel barriers. Also, these authors did not provide 
a symmetrized formula in the steady state, which is essen- 
tial for the calculation and discussion of the spin current. 
Mu et ali^ used the difference between the charge cur- 
rents through the spin up and down channels to define 
the z-component of the spin current. However, these au- 
thors did not properly take into account the difference of 
the two local quantization axes of the two ferromagnetic 
leads which strongly affects the tunneling hamiltonian as 
pointed out by Rudzihski et al."^^ Moreover, their expres- 
sion of the charge current was not symmetrized. As a 
result, their formula is correct only for the parallel case. 



III. CALCULATION OF THE KELDYSH 
GREEN'S FUNCTIONS 



To investigate the nonequilibrium transport proper- 
ties, there are two commonly used techniques to calcu- 
late the Keldysh Green's functions. One is the real-time 
diagrammatic techniquef^ '^^i^^'^^i'^^ based on a perturba- 
tion expansion in terms of the dot-lead coupling strength, 
whereas the Coulomb interactions on the dot are exactly 
taken into account. However, this technique only con- 
siders finite order tunneling processes, and cannot deal 
with the coupling between the dot and the electrode ex- 
actly. The other alternative is the equation-of-motion 
technique^iiSii2ii^i2ii22 which treats the dot-lead coupling 
exactly, while the strong correlations on the dot can be 
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dealt with only approximately. 

In this paper, the Green's functions are solved by the 
equation-of-motion technique with the Hartree-Fock de- 
coupling schemei^SiSl The solution can be written in a 
compact form of the matrix Dyson equation 



Gh iuj) 



where 



gd (w) 



{<^-(^d)i'^-<^d-U) (u)-ed){uJ-<id-U) 
UjnT.i) ^-ed-[/(l-(n-r.T)) 
(w-Cd)(w~ed-;7) (uj-£d)(u)-ea-U) 



(12) 



with {ria^p) — {d]^df)) and gd (w) denotes the correspond- 
ing Green functions in the matrix form of the uncoupled 
dot. The self-energy S^'^^ [uj) is given by 



with 



.(0) 



^ - ek,±,L 

2 



K 



k,i? 



I?K sin ( 



K 



cos2 (61/2) sin^ (61/2) 

^ — ek,±,i? ^ — ek,^,i? ' 
1 1 



(13) 



w — ek,+,-R ^ ^ £k,-,_R 
Then one can calculate the retarded Green functions as 

.(0)r , 



(C.) 

Gli (^) 



{lo) - A(c.)SL^i' (c.)J lB{u), 
'gl, (.;)-yl(c.)sVtHl /SH, 



where 

A{lj) = g'^^ (uj) gl^ {uj) - g\^ (w) gl^ {uj) , 

B{oj) = i-.9[^(^)sVtH-.9liMsL"^''H 



gl^ (^)sLt(c^)-5lT MSr'M 



(0)r 



+A(a;) 



.(0)r , 



(w) SL"!" (co) - E^'r (co) E^Y' C'^) 



(O)i- 



.(0)-r 



The retarded self-energies ^±± (w) and S^%"(a;) are 
given by the formulas 

sitM - -^[r±,i(c.) + r±,«(.^)cos2(0/2) 



.(0)r 



(CO) 



-fr^,fl,(^)sin2 (6/2)], 

1 

-7 [r+,i?(t^) - r_,_R(tj)] sin^ 



(14) 



In the following we assume 

r±,L(t^) = r±,L = To (1 ±pi) , 
r±,fl(t^) = r±,_R = (1 ±pr) , 

where and pr denote the spin polarization of the left 
and right electrodes, respectively, and the parameter 7 
expresses the asymmetry coupling of the left and right 
electrodes to the dot. {uj) can be obtained from the 
Keldysh equation. 



G< (c.) = G^ (c.)S<(c.)G5 (co), 



(15) 



where the full self-energy S< {uj) is related to S(°)< (uj) 
via the Ng ansatz'^* 

^^±{UJ) = lTo[fL{uj)il±Pl)+7fR{uj)il±PrCOS0)], 

^±T (^) = ijTofji{uj)prSm6. (16) 

The statistical averages of (ua^p) have to be calculated 
self-consistently in the following way: 



(ric^a) Im y 



duj 
2^ 
duj 
2^ 



(uj) 



(uj) . 



(17) 



This approximate calculation of the Keldysh Green's 
functions does not take into account the Kondo-like 
correlations which need a careful treatment of the 
Coulomb interaction on the dot. Some previous 
works which studied the charge transport properties 
of this system discussed the Kondo effect, including 
the collinear alignment^'^'^'^ ' lOiHii^iis^ie g^^^^ ^j^g non- 
collinear case^i^iiS^iSl. It is left for our future work to 
discuss the influence of the Kondo-like correlations on 
the spin transport properties of this system. 



IV. RESULTS AND DISCUSSIONS 

Now we numerically calculate the three components of 
the spin current. Since a general magnetic configuration 
of the leads is considered, the spin tunneling magnetore- 
sistance (STMR) can be estimated by 

STMRa^ — ^ ,a = x,y,z (18) 

where Jx,y,z{&) denote the three components of the spin 
current. In the following three different situations are 
considered: a symmetric junction with fully polarized 
external electrodes {pi = Pr — 1), with partially polarized 
external electrodes {pi = p,. = 0.4), and an asymmetric 
junction {pi — 0.4, = !)■ 

The J^-voltage curve for the symmetric cases reveals 
typical step-like characteristics. Below the lower thresh- 
old voltage, the dot is empty and the sequential contri- 
bution to Jz is exponentially suppressed. The first step 
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in Jz occurs at a critical bias, where the discrete level 
ed crosses the Fermi level, whereas the step at a higher 
threshold corresponds to the case when + C/ crosses the 
Fermi level. In the same voltage range, Jz in the case of 
Pi = = 1 is much larger, since the external electrodes 
arc fully polarized. Mu et al?~ also considered this case 
{pi = Pr = 0.4,6 = 7r/3), but their result is different 
from ours because they did not properly take into ac- 
count the difference of the local quantization axes in the 
two leads and hence the result for the ^-component of 
the spin current is incorrect. The case of pi = 0.4, = 1 
is more complex (red dashed curve in Fig. 1(a)), as the 
asymmetry between the left and right electrodes gives rise 
to asymmetrical transport characteristics of the junction 
with respect to the bias reversal. For the positive bias, Jz 
curve is rather smooth above the first threshold voltage, 
while for the negative bias, below the first threshold se- 
quential tunneling is exponentially suppressed and only 
the higher-order tunneling processes are possible. When 
Cd approaches the Fermi level of the left electrode, the 
resonant tunneling can happen. However, as the bias 
further increases, Jz is suppressed by an electron on the 
QD since the electrode is partially polarized (Coulomb 
blockade effect), and a small peak appears as a result 
of competition between the resonance tunneling and the 
Coulomb repulsion. After the second resonant tunnel- 
ing, Jz finally saturates at a certain level. The behavior 
of spin current component J^ (Fig. 2(b)) is similar to the 
component Jz (Fig. 2 (a)), because the magnetizations of 
the two leads are aligned in the x-z plane. However, the 
asymmetry effect resulting in the appearance of a peak 
at the first threshold is more pronounced. It appears 
even for the symmetric electrodes {pi ~ Pr — 0.4, 1), be- 
cause the Coulomb blockade effect already shows up. The 
asymmetry of the spin current curve is even more pro- 
nounced for the y— component (Fig. 2(c)). Nevertheless, 
the two peaks on the Jy curve are exactly located at the 
two resonant tunneling biases. 

The spin current is strongly affected by the angle 9 be- 
tween the magnetic moments of the leads and we can 
use STMR to describe it. In the free regime, where 
\eV/2\ > Ed + U, the QD energy level may be occupied 
by two electrons, because the Coulomb correlation plays 
a little role in the spin tunneling. As a result, Jz and 
STMRz exhibit a monotonic variation between the par- 
allel and antiparallel magnetization configurations, which 
is typical of a normal spin- valve effect. Under the third 
condition [pi = 0.4, = 1), Jz can achieve a negative 
value. The absolute values of Jx and Jy achieve their 
maxima between 6 — and 6 = n as shown in Fig. 3(c) 
and Fig. 3(e), since the absolute values of the x and y 
components of the electron spin in the right electrode 
may increase when the magnetic moment of the right 
lead approaches the x-y plane. 

In the Coulomb blockade regime < \eV/2\ < ed + 
U, the QD energy level can be occupied only by one 
electron. The Coulomb interaction plays an important 
role in the spin current through the QD. In Fig. 4(a), it 
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FIG. 2: Voltage bias dependence of the spin current for 
9 = 7r/3. (a) Jz, (b) J^c, and (c) Jy. The parameter values 
assumed are: ea = O.leV, U = 0.4eV, To = O.OleV, 7 = 1, 
and r = lOOA'. 



is found that Jz (6 = 0) is no longer maximal and Jz (9) 
is greater than Jz {6 = 0) in a wide range of 9 under this 
condition {pi = O.A,pr = 1). It is quite different from that 
in the free regime. The coupling between the QD and 
the ferromagnetic leads may induce an effective exchange 
field, and its strength and orientation with respect to the 
global quantization axis depend on the bias voltage and 
the angle between magnetizations of the leads. When 
only one electron resides on the QD energy level, the 
spin degrees of freedom experience a torque due to the 
effective exchange field, which results in precession of the 
spin around the field. This process would suppress Jz, 
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and the competition of the spin precession effect and the 
spin- valve effect leads to the anomaly of Jz {9). As a 
result of the spin precession, the signs of Jx (6) and Jy {6) 
are opposite to those in the free regime. 

In conclusion, we have derived a general formula for the 
spin current through a QD coupled to ferromagnetic leads 
with noncoUinear magnetizations, and used the formula 
to calculate the spin transport properties of the system. 
The competition of the spin precession and the spin- valve 
effect results in an anomaly of the angle dependence of 
the spin current. Further investigations are needed to 
carefully treat the Coulomb interaction on the QD. 
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